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In a recent letter Kumar et al. outline a model for thermal transport in nanofluids that purports to explain 
the anomalous thermal conductivity enhancements observed experimentally. Since this effect is expected to have 
significant technological consequences and may already be relevant for certain applications 2] this is an important 
exercise. Unfortunately, the model of Kumar et al. falls short on several conceptual counts. 

The authors start by considering the case of fixed solid particles and postulate the existence of "two parallel paths 
of heat flow" , one through the suspending liquid and the other through the dispersed solid. This is a questionable 
assumption, particularly for the extremely dilute suspensions against which they wish to test their model. By pursuing 
this argument along an approximate derivation Kumar et al. produce a relation for the effective thermal conductivity 
of the nanofluid, Eq. 7 of |f(, which includes a dependence on the radiuses of the solid particles, r p , and liquid 
molecules, r m : 
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,where k p and k m are the thermal conductivities of the solid and liquid respectively, and e is the volume fraction of 
the dispersed solid particles. One problem with the above relation is immediately apparent if we notice that at fixed 
e, r m /r p — > implies k e ff/k m — > 1, i.e. the solid inclusions have no effect on the thermal conductivity if they are 
much larger than the liquid molecules. In fact, this situation is rigorously described by theories that treat the liquid 
as a continuum and where accurate, widely accepted results are available for all e 0. The dependence of Eq. I on 
kp/km is very different from these results even at r m /r p > 0. For example, if the solid clusters are much less thermally 
conducting than the liquid, k p /k m — > 0, Eq. I predicts k e ff/k m — > I, while if they arc much more conducting, 
kp/k m — > oo, it yields k e ff/k m — > oo, both independently of e (and r m /r p ), and both of which can only be deemed 
incorrect. 

A dependence of & e // on r p can arise for example due to liquid-solid interface thermal resistance Q, even when 
the liquid can be treated as a continuum. Furthermore, as r p /r m decreases the effective volume fraction of the solid 
particles, e.g. approximated by e e ff = e(l + r m /r p ) 3 for e <C 1, should also perhaps replace e in conductivity estimates, 
introducing an additional dependence on r p . However, Eq. I captures no such effects and fails some simple tests. 

Kumar et al. consider next the effect of the solid particles motion. To this end they invoke kinetic theory and 
introduce cu p as the "thermal conductivity of the particle" , where u p is an "average particle velocity" . Their complete 
thermal conductivity model is obtained by replacing k p in Eq. I with cu p . This is quite problematic since it has 
the effect of eliminating k p as a parameter in determining the nanofluid thermal conductivity, and thus makes the 
previous analysis for fixed solid clusters rather futile. For example, if u p — » 0, .e.g. the suspending liquid is frozen, 
k e ff — > k m , i.e. the solid inclusions have no effect on fc e // irrespective of k p and e, when this case should reduce to 
the fixed particles problem. 

In fact, cu p cannot be interpreted at all as the thermal conductivity of a solid particle. By the authors own kinetic 
theory arguments ki = cu p is the thermal conductivity of a dilute gas of particles that possess internal energy 
In principle, a quantity like ki could provide an estimate for direct, Brownian motion transport of heat by the solid 
clusters, but it should not supplant k pi the thermal conductivity of the solid. Unfortunately, the ki calculated in 
Q is inadequate even for this purpose. It is not clear why the authors definition of u p would be appropriate in the 
given context, while their estimate of the "constant" c oc nlc v (which is not dimensionless as conveyed in the paper), 
with n - number density of the solid particles, I - their "mean-free path" and c v - heat capacity of a solid particle, 
apparently assumes I cx l/nd p , d p = 2r p , which only takes into account the rare collisions between the solid particles 
and ignores the effect of the liquid. A more reasonable estimate for ki is ki ~ nDc v , with D = ksT '/3Trr)d p - Stokes- 
Einstein relation. This yields for d p = lOnm gold particles in water at normal conditions and e = 1%, ki/k rn ~ 10 -6 . 
Therefore, the direct Brownian motion contribution to thermal transport can be safely ignored for nanofluids. Finally, 
the theoretical values plotted by the authors in their Fig. 4 are off by orders of magnitude from their own formula 
Eq. 10, purportedly used to calculate them. 

This work was performed under the auspices of the U. S. Department of Energy by University of California Lawrence 
Livermore National Laboratory under Contract No. W-7405-Eng-48. 
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